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Abstract —The Teager-Kaiser energy operator (TKO) belongs 
to a class of antocorrelators and their linear combination that can 
track the instantaneous energy of a nonstationary sinusoidal sig¬ 
nal source. TKO-based monocomponent AM-FM demodulation 
algorithms work under the basic assumption that the operator 
outputs are always positive. In the absence of noise, this is assured 
for pure sinusoidal inputs and the instantaneous property is also 
guaranteed. Noise invalidates both of these, particularly under 
small signal conditions. Post-detection filtering and thresholding 
are of use to reestablish these at the cost of some time to 
acquire. Key questions are: (a) how many samples must one 
use and (b) how much noise power at the detector input can 
one tolerate. Results of study of the role of delay and the limits 
imposed by additive Gaussian noise are presented along with the 
computation of the cumulants and probability density functions 
of the individual quadratic forms and their ratios. 

Keywords — Teager-Kaiser energy operators; autocorrelators; ra¬ 
tios of quadratic forms 


I. Introduction 

There are many applications in communication and control 
where decisions are made on the basis of the ratios of quanti¬ 
ties, measured or computed, which are compared with thresh¬ 
olds. The present work is concerned with situations, where 
variables are outputs of a linear combination of quadratic 
forms in correlated Gaussian random variables and their ra¬ 
tios. Such quantities arise in a variety of contexts including, 
for example, monocomponent AM-FM demodulation using a 
special class of quadratic detectors called the Teager-Kaiser 
Energy Operator (TKO), maximum ratio diversity combining 
schemes to combat fading, directivity optimization of antenna 
arrays, calculations of SINR and various error and outage 
probabilities, etc. This work primarily examines the noise 
sensitivity of the TKO by studying the statistical distribution 
of the operator outputs and their ratios. 

The TKO is an important class of product correlators used 
for nonstationary signal analysis. Starting in the early 1980s, 
Herbert M. Teager and Shushan M. Teager [1] developed 
several tools, including the energy operator in a series of 
experiments for modeling nonlinear structures in speech sig¬ 
nals, caused by modulations and turbulence. Following their 
pioneering work, the operator was systematically introduced 
by Kaiser [2] and was shown to track the energy of a simple 
harmonic oscillator. The past two decades have witnessed 
extensive research showing the usefulness of the operators in 
diverse fields, including AM-FM decomposition for monocom¬ 
ponent signals [3], [4] and aided by the use of filter banks for 
multicomponent signals [5], speech analysis [4], image texture 
analysis [8], cochannel and adjacent-channel signal separation 
[7], frequency and time discrimination [9], [10], etc. to name 
a few. 


For continuous-time signals, x{t), y{t), the kernel [x,y] = 
xy — xy measures the instantaneous differences in the relative 
rate of change between x and y [11], where the overdots repre¬ 
sent time derivatives. If x is the Hilbert transform of y so that x 
and y form an analytic pair, [x,y] gives the intensity-weighted 
phase derivative and the instantaneous frequency (IF) is defined 
as uj{t) = ■ if y = X, then we get the continuous-time 

TKO: = [x,x] = x^ — xx. For a monocomponent non¬ 

stationary sinusoid x{t) with time-varying amplitude a{t) and 
instantaneous frequency the operator applied to the signal 
and its derivative yield respectively, 'k[a:(f)] « a'^{t)uj^{f) and 
~ a?{f)uj'^{f), where the approximations hold if the 
bandwidths of the amplitude/frequency modulating signals are 
reasonably small compared to the carrier frequency. Under the 
assumption that the operator outputs are always positive, these 
observations have motivated the continuous energy separation 
algorithm (ESA) [3], [4] for AM-FM decomposition. 
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The discrete-time version of the operator is derived by 
replacing x{t) with sequences x[n] = x(nT) and derivatives 
X with two-sample backward differences (x\n] — x\n — f])/T 
where T is the sampling interval. The continuous-time operator 
then reduces (up to one sample shift) to the discrete version 
given by T'[a:[n]] = {x^[n] — x[n — l\x[n + 1])/T^ [2]. We 
introduce a general four-sample symmetric kernel as a linear 
combination of two correlators, viz., = {x[n—p]x[n+ 

p]—x[n — q]x[n + q])/T'^. Forp = 0, g = 1, T*® reduces to the 
conventional discrete-time TKO T'q [2]. The effect of applying 
T'g at a sampling rate of fs = 1/T is equivalent to applying thj 
at a sampling rate of fs/q, so that q defines the underlying sub¬ 
rate processing. Replacing the forward and backward sample 
delays with time delays, and ignoring the scaling factor T, 'F® 
as a function of the delay parameters {tp,tq) becomes 


'f^lix) = x{t - tp)x{t A tp) - x{t - tq)x{t + tq), p < q (2) 

A distinct advantage of the TKO is its superior localization 
property. When the input to the operator is the sum of a deter¬ 
ministic signal and a Gaussian noise process with a specified 
covariance function, the discrete operator kernel is an indefinite 
quadratic form in nonsingular normal vectors. In the absence 
of noise, conditions for the positivity of the operator output can 
be expressed in terms of the concavity of the logarithm of the 
signal magnitude [6] that hold in general for pure sinusoids and 
the instantaneous property is also guaranteed. Noise invalidates 
both of these, particularly under small signal conditions. The 
frequency computation involves noisy differentiation opera¬ 
tions and a short time window that precludes any scope for 





averaging. Because of the inherent nonlinearity, significant 
noise X noise and signal x noise perturbations exist. The effect 
of these terms becomes pronounced when squaring or ratio 
operations are employed to estimate signal parameters, and a 
suitable post-filtering action is necessary at the cost of some 
time to acquire. For the existence of moments of the ratio, the 
behavior of the density of the denominator near singularities 
plays a vital role. Under small signal conditions, large am¬ 
plitudes can occur frequently when the denominator assumes 
infinitesimally small values, so that the integral defining the 
moments diverges because of long tails in the density function. 
A way out is to eliminate negative outputs and outputs near 
zero by conditioning all measurements of the ratios on the 
event that the denominator exceeds some detection threshold— 
the measurement is ignored if threshold detection fails to occur, 
thus precluding a division by zero [12]. Thresholding retains 
the speed advantage of the operator and is a viable alternative, 
when filtering is not an option. 

Earlier work has examined the problem of noise and found 
expressions for the variance of the amplitude and frequency 
estimates. Kaiser [2] briefly mentioned results for the discrete 
operator in zero-mean additive white Gaussian noise. Bovik 
and Maragos et al. [5] gave a more extensive treatment of 
multiband demodulation and developed the basic statistical 
properties of the TKO. A spline-based ESA was proposed in 
[13] for robust estimation in the presence of noise. A more 
recent work on short-term energy estimation in additive noise 
is given in [14]. Statistical distributions of the operator outputs 
and their ratios have not received any attention. This work is 
primarily aimed at filling this gap by evaluating the cumulants 
and probability density functions of the individual operator 
outputs and their ratios. 

II. Properties of the TKO in Relation to 
Erequency Discrimination 

The symmetric discrete TKO operator kernel implements 
a translation-invariant, sliding-window filter that integrates 
energy detection with frequency filtering [15]. The frequency 
responses of the three-sample input operator and four- 
sample input operator T'* have an important difference. Eor 
a single sine wave of unit amplitude and frequency w, the 
frequency response 4'q(u;) = 1 — cos(2wfg) is unipolar 
and periodic. At low frequencies (wT < 7r/4q) when the 
linear approximation sin(u;T) « wT holds, oc uP. 

Traditionally, this range is used for AM-EM demodulation. 
Eor a periodic input signal, it is easy to verify that the power 
in the odd harmonics can be obtained by averaging the output 
of a three-sample input TKO T'g by adjusting the delay tq to 
correspond to a phase difference of 7r/2 for the fundamental 
frequency. 

On the other hand, the frequency response = 

cos{2(jjtp) — cos{2ujtq) = 2s\nuj{tq — tp){tq + tp) is bipolar and 
is not proportional to at low frequencies. The response is 
zero at w = 0 and when either oj{tq — tp) or uj{tq+tp) is a mul¬ 
tiple of TT. One observes regions where the response is linear 
around the zeros of The location of the extrema of the 

frequency response can be identified as a solution of the equa¬ 
tion, tpSm{2ujtp) = tqSm{2ujtq). A discriminator characteristic 
thus results. When the input consists of several sine waves of 
amplitudes ak and frequencies ujk, i.e., x{t) = ^afcsin(u;fcf), 
one gets a dc term, ^a|[sin^(wfcfq) — sin^(wfcfp)] and a time- 
varying cross-term, which is a combination of the sum and 


difference frequency terms arising out of mixing operations. 
A typical cross-term component is of the form (for k < ly. 
akai{cos{uJk-uJi)t[cos{uJk+uJi)tp-cos{uJk+uJi)tq]-cos{uJk + 
uJi)t[cos{uJk — 0 Ji)tp — cos{(jJk — uji)tq]}. The choice of a 
symmetric kernel ensures that all second harmonic compo¬ 
nents get cancelled. When the input is differentiated to give 
x{t) = ^a;fcafcCOs(wfcf), all the cross-term components are 
multiplied by WfcW;. Ordinarily, for monocomponent signals, 
the dc term obtained from the TKO for the signal and its 
derivative are used to find estimates of IE and amplitude. Eor 
demodulating multicomponent signals, filter banks are used 
[5]. When two frequencies are present simultaneously, the 
T'® outputs at the sum and difference frequency cross-terms 
are of equal amplitude. If the ratio of these frequencies is 
determined, one can utilize this information in combination 
with the and 4'9(i;) outputs to estimate amplitude and 

frequency parameters. So, if the signals are closely-spaced [9] 
or belong to the same filter bank, frequency discrimination 
using the sum and difference frequency cross-terms generated 
by a suitable combination of correlators is a viable alternative. 

III. Experiments: TKO Outputs for a Sum of Two 
Sinusoids and Sinusoid Plus Noise Input 

Given the wide usage of TKOs in IE estimation for 
monocomponent signals, it is important to assess the effect 
of noise or of an interfering signal upon the IE estimates, and 
check for any physical inconsistencies thereof. The situation is 
best analyzed by considering the effect of a sinusoidal inter¬ 
ference. Considering the most general form of a sum of two 
sinusoids with positive amplitudes oi and 02 , incommensurate 
frequencies fi and / 2 , and independently uniformly distributed 
phases cpQ and yjo, 

x(f;ai,/i, 025 / 2 ) = aiC0s(27r/if -f ipo) + a2C0s(27r/2f -f i/o) 

= aiCOS(/? -f a2C0St/> 

Without loss of generality, the model can be simplified by 
taking into account the amplitude and frequency ratios, a = 
02 / 01 , / = / 2//1 and Oq = ipQ — (fiQ as the following analysis 
is sensitive only to such relative parameters. The simplified 
signal model becomes 

x(t;a,f) = cos(27rf) -f acos(27r/f -f Oq) (3) 

Eor certain parameter ranges in the a-f plane, the 
continuous-time TKO 'I'(x) = x^ — xx can go negative, 
a situation that clearly undermines the physical significance 
of the operator as an energy detector. The situation is best 
analyzed at the extremas of x{t), where the probability of such 
an event occurring is the highest. At an extrema, the derivative 
of (3) is zero, i.e. x{t;a,f) oc sin(27rf)-|-a/sin(27r/f-|-0o) = 0. 
The conditions for which the operator output is negative at such 
an extrema is found by evaluating the condition, —xx < 0, i.e., 

a^/^cos^(27r/f -f Oq) + cos^(27rf) 

+ l)cos(27rf)cos(27r/f -f 0o) ^ 0 (4) 

The required condition is obtained by solving the inequality (4) 
which is a quadratic in cos(27r/f + Oq), giving the following 
result: If x{t;a,f) admits an extremum at t = to, whenever 
cos{2Trft + Oq) evaluated at such an extrema lies between the 
curves yR = —cos(27rf)/a and yc = —cos(27rf)/a/^, T'(a;) 
will be negative. Of course, the operator output can be negative 
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Fig. 1. TKO Response for two sinusoid input: a = 0.6,/ = 2.3. x{t) = cos(27ri) + acos(27r/f + 9o). In the top panel, the derivative of the first component 
(Ml) and the opposite of the derivative of the second component (M 2 ) are plotted, so that their intersections correspond to the extremas of x(t). The marked 
extrema locations contribute to the zero-crossings of the TKO output shown in the bottom panel. The middle panel shows that whenever cos(27r/f -F 0o) 
evaluated at the extrema points lies between the curves marked yQ and yn, the corresponding extrema of x{t) would contribute to the zero crossing of the 
TKO output 


at points other than the extrema, albeit with a lesser probability. 
The situation is depicted in Fig. 1. 

The foregoing analysis for the sum of two sinusoids is 
easily extended to the more general, sinusoid plus noise input 
case, by considering a spectral distribution of the noise concen¬ 
trated solely at a frequency /, making the noise a sinusoidal 
Gaussian process with Rayleigh-distributed amplitude [16]. 
Intuitively, the most vulnerable are the inflection points given 
by the positive minimas (x > 0, i = 0, x > 0) and the 
negative maximas (x < 0, i = 0, x < 0). The zero-crossing 
rate of the operator output as well as the distribution of the 
time duration for which the operator output remains negative 
provide valuable insights into the problem at hand. In partic¬ 
ular, the average time duration for which the operator output 
remains negative can provide valuable insight in deciding on 
appropriate delay parameters (e.g., tp, tq for T'P or a suitable 
integration time constant in the post-detection filtering stage. 

IV. Statistical Distributions of TKO Outputs 
AND Their Ratios 

A. Probability Density Functions of the TKO outputs 

Notation: tr(-), (•)', (•)^ denote respectively, the trace, 
nonconjugate and conjugate transpose operators. E[-] is the 
expectation operator and Im[-] is the imaginary part operator. 
i = is the imaginary unit. The characteristic function 

(chf) of a random variable (RV) X is 4>x(f,) = E[exp(i^X)]. 
We write fx{x) = Pr{X = x} to denote the probability 
density function (pdf) and Fx{x) to denote the cumulative 
distribution function (cdf) of RV X. 


We first consider the case when the signals at the detector 
input are contaminated with additive, zero-mean, stationary 
Gaussian noise with power spectral density Nq. Let X be the 
column matrix formed from r real signal plus noise samples 
Xj, j = Let pL = EX be the matrix of the means and 

M = E[(X —EX)(X —EX)'] the positive definite (unnormal¬ 
ized) covariance matrix, i.e., X ^ Afr{pi,M). For the jointly 
Gaussian random variates (x[n — q\,x[n\,x[n + <7]), the matrix 

/ 0 0 -i\ 

of the symmetric kernel T'l] is given by J = I 0 1 0 ). 

v-5 0 0 y 

Likewise, for jointly distributed RVs (x[n —q\^x[n —p\,x[n + 
p],x[n + gj), the matrix of the symmetric kernel T'* is given 

( 000 

olio]- 

-i 0 0 0 / 

Equivalently, the TKO outputs can be viewed as a linear 
combination of quadratic detectors with the product correla¬ 
tions reduced to a sum and difference of squares. Thus, for 
instance, for jointly distributed Gaussian RVs (x[n],s[n],(i[n]) 
where s[n] = (x[n + g] + x[n — q])/2 and d[n] = (x[n + 
q] — x[n + g])/2, = x'^[n\ — s^[n] + (P[n] and 

J = diag(l,-l,l). 

Given J, the related quadratic form in X is V{X) = 
X'JX. For a zero-mean, stationary random process whose 
spectrum is Gaussian, the covariance kernel may take the form 
R{t) = exp(— cf^). Letting M = LL', Y = L~^X and 
Z = Y — L~^p,, we have Z ~ Nr{0,I), where I is the 
identity matrix. Let P be a r x r orthogonal matrix which 



























diagonalizes L'JL, so that P'L'JLP = diag{Xi,...,Xr) = 
A and PP' = I, where Ai,.--,Ar are the eigenvalues of 
MJ. Then, letting U = P'Z, we have U ~ Afr{0,I) 
and V{X) = {Z + L-^g)'L'JL(Z + = (U + 

s)'k{U + s) = Y.]^iXj{Uj+Sjf = where 

s = P'L~^g, = {si,...,Sry and Wj s are independent RVs 
with a single degree of freedom and noncentrality parameter 

s|, i.e., ~ For an indefinite quadratic form where 

the A's can be any real number, letting Xj > 0 for j = 

Xj < 0 for j = k + and Xj = 0 for j = / + it is 

easy to see that V = Vi — V 2 , where Vi = '^j^iXjWj and 

F 2 = - Thus, when X ~ KipM), M > 0, 

indefinite quadratic forms in X are distributed as the difference 
of two weighted sums of independent x^ RVs [17]. 

The characteristic function (chf) of V{X) is given by 




exp{-i/r'[7 -{I- 


(5) 


Alternatively, (5) can be put into the following diagonal form 
[27], [17], [18], 


= 11^^“ 2i^XkNo) ^^^exp{ 
i=i 


s'j 2i^XjNo 
2No 1 - 2i^XjNo ^ 


(6) 


For the more specialized case of narrowband inputs, the 
correlator outputs can be represented as a linear combination of 
independent, scaled non-central x^ variables with two degrees 
of freedom [18], [19]. Now the noise power per cycle Nq is 
equally split between the in-phase and quadrature components, 
so that for each in-phase component of noise at the detector 
input, an additional independent quadrature component of 
noise is added to the output. The latter contributes an extra 
signal independent term so that in effect, the chf in ( 6 ) is 
multiplied by its value for noise alone and we get [26], [27], 

= n(l - ^CAfcA^o)~^exp{^ ^ 

Ao 1 - i^XjNo 


Formally, let C be the column matrix formed from r complex 
RVs, Cj = Qj +ibj, j = l,...,r, where aj, bj are normally dis¬ 
tributed RVs with E[(afc—Eafc)(ai—Ea;)] = E[{bk — 'Ebk){bi — 
E 6 i)], and E[(afc-Eafc)(&i-E 6 ;)] = -E[{ai-Eai){bk-Ebk)] 
Let EC = C and L = E[(C-EC)(C-EC)t] be the complex 
(nonsingular) covariance matrix. Then, given any Hermitian 
matrix Q, the chf of the (real) quadratic form V{C) = C^QC 
is given by [19] 




exp{-Ct[/ -{I- 

\I-t^LQ\ 


( 8 ) 


when a straightforward calculation of the matrix operations 
leads to (7). 

In general, inversion of the chf does not yield closed-form 
solutions for the pdf, except for a few trivial cases. For a single 
eigenvalue A, using the Laplace transform pair 29.3.77 [20], 
inversion of the chf in ( 6 ) yields [27] 


fviv) = 


exp{-2^[g" + x]} 

\/2'kNqXv 


cosh 




, v > 0, 


and zero otherwise. Likewise, for inversion of the chf in (7), 
the equivalent expression is that of a Rician distribution for 
the envelope [26]. 


fviv) = 


exp< - 


{ No + a]} 


XNn 




and zero otherwise, where Iq{-) is the zeroth-order modified 
Bessel function of the first kind. 

For positive only eigenvalues, one can use non-causal 
convolution of the individual density functions corresponding 
to the independent eigenmodes yielding the composite density 
function /vx(u+), and similarly for negative only eigenvalues, 
the density function, /y_(u_). The final density function 
fviv) of V is obtained by convolution of the independent 
densities. 


00 

fviv)= J dv+fv+iv+)fv.iv-v+) (9) 

max( 0 , 7 i) 


If V > 0 , the lower limit of is v, while if v < 0, the lower 
limit is zero. 

As an example, consider the central case (/r = 0) with 
eigenvalues Ai, A 2 >0 and A 3 , A 4 < 0. The corresponding 
chfs are given by, 

cj)v+ = (1 - 2j^XiNo)-^/H 1 - 2j^X2No)-^/^ 


(t>v- = (1 + 2jeA3iVo)-^/2(l + 2jeA4iVo)-'/2 


Using the Laplace transform pair 29.3.49 [20], the correspond¬ 
ing density functions are 


fv+iv+) = 




for > 0 and zero otherwise. 


fv-iv-) 


for u_ < 0 and zero otherwise. Jo( ) is the zeroth-order Bessel 
function of the first kind. The final density function can be 
found by convolving /vx(u+) and /v_(v_) using (9). 

Efficient numerical algorithms exist for evaluating the cdf 
of a quadratic form [21], which mostly rely on Gil-Pelaez’s 
inversion theorem [ 22 ]: 


0 

The cumulant generating function admits a power series 
expansion KyiC) ~ log(/)y(^) = cumu- 

lants Ks are defined by the following identity [23]: 

Ks = 2«-^(s - l)!{tr(MJ)« -f sg'JiMjy-^g} (11) 


The cumulants ki and K 2 are the mean and the variance, 
respectively. For a pair of independent RVs (X,V), KgfpX + 
qY) = p^KsiX) -f q^KsiY). For s > 2, the standardized 
cumulants given by ps = are invariant under affine 

transformations and provide summary measures of departure 
from normality in the sense that for a normal distribution, 
Ks = 0 for s > 2. In particular, p^ measures the skewness 
and p 4 measures the kurtosis of the distribution. 



B. Probability Density Functions of the Ratio of Two TKO 
Outputs 

The statistical properties of ratios of quadratic forms, 
their moments, densities and distributions, or approximations 
thereof, have been extensively studied [23]-[25], [29]. The 
reader is referred to [17] for a survey of the mainstream 
statistical literature. In general, for any two real functions Vi 
and V 2 (which need not be quadratic forms), the distribution 
of the ratio, R — ^ can be computed if the joint density 
function of Vi and V 2 , /vi.vb is known. Assuming that 

V 2 is not negative, the probability density of R is given by 

00 

fR{r)= Jv2fVuV2i'^'^2,V2)dV2. 

0 

However, it is very rarely that the joint density is known or Vi 
and V 2 are independent in which case a simple relation applies. 
When Vi and V 2 are not necessarily independent, Geary [28] 
showed that the density of the ratio R can be written as 


fnir) 


00 



— 00 


d^2 


dil 

i2 = -rii 


( 12 ) 


When Vi and V 2 are noncentral quadratic forms in cor¬ 
related Gaussian RVs, an efficient numerical algorithm is 
given in [29] that uses Geary’s formulation for computing the 
density function of the ratio R. The technique can be used for 
computing the pdf of the ratios and for the 

instantaneous frequency squared and instantaneous amplitude 
or envelope squared estimates, respectively. For the latter, 
the density function for the numerator which is of the form, 
V 3 = Vi can be obtained from the density function of Vi 
by a nonlinear transformation. Correspondingly, /vj 14 (^ 3 )^ 2 ) 
needs be transformed for use in ( 12 ), to get the pdf of the ratio 

V^IV2. 


V. Results and Discussion 
An important consequence of the assumption of a sta¬ 
tionary Gaussian noise structure at the detector input is that 
all higher order probability densities are completely speci¬ 
fied if the covariance function of the input noise is known. 
Because of stationarity, the covariance matrix has elements 
rriij = 'KIxiXj], which depend only on the absolute difference 
or distance |z—jj between the samples and accordingly, a single 
subscript rrij = WyXiXj], i,j = 0,±1,±2,... etc., suffices, where 
Xi = x[i]. The covariance function of the input noise process 
is taken to be a Gaussian kernel given by exp(—f^/2). Smooth 
estimates of the signal derivative are obtained from a finer grid 
by interpolating the signal samples. The first few moments of 
the TKO output are easy to find. For instance, writing the 
output of vk® as (si -I- ni)(s2 + 112) - (sa + n3)(s4 + nfj, 
where Si and rti are the signal and noise components at sample 
point i, we have ]E[T'9] = (siS 2 - S 3 S 4 ) + {m 2 p - m 2 q)- 
The variance Var['I'®] has the form + 2 [siS2 TO2 p+ 

S3S4TO2q - (S 1 S 3 + S 2 sf)mq-p - (siS 4 -f S 2 S 3 )mq+p\. We 
see that there is an average term contributed solely by noise, 
as well as a few signal x noise terms that cause strong time- 
varying fluctuations. For a time-varying signal, the mean, 
variance, and the covariance are also time-varying and an 
averaging over one complete cycle of the input sinusoid 
is necessary to have meaningful expressions for the mean. 


X 10 ^ Teager Output Pdf for signal: T(x) 



X 10 ^ Teager Output Pdf for signal derivative ; 'T(:ii:) 



TKO Output (normalized by input noise power) 



Ratio R 

Fig. 2. Pdfs for the TKO outputs (normalized by the noise power into 
detector) for the signal, signal derivative and the ratios (Input SNR =17 dB). 
For the distribution of the ratios, two different sets of pdfs show the effect 
of noise power into detector: the sets (a,b,c) and (a’,b’,c’) correspond to an 
input SNR of 17 dB and 25 dB respectively 


variance as well as the higher-order moments. Also, signal and 
signal derivative product terms lead to the creation of second 
harmonics, and a suitable post-filtering action is necessary. 

The computed probability densities of the outputs of dif¬ 
ferent TKO variants, as well as the ratio which yields 

the instantaneous frequency squared are shown in Fig. 2 for 
different input signal-to-noise ratios (SNRs). As expected, the 
mean output SNR IE[Ki(s_|_iv)/fJi(Ar)] improved as the delay 



















































Fig. 3. Pdfs for the individual quadratic forms: (/ 1 Q 2 — l 2 Ql), {I 1 Q 2 + 
I 2 Q 1 ), and (7j + + Qf + Q^)- Input SNR = 9.04 dB. Inset: pdfs of the 

sine and tangent I-Q correlators (a) 2{IiQ2 — I 2 Q 1 )/ +1^1 + Ql + Q 2 ) 
and (b) (I 1 Q 2 — hQl)/(I 1 Q 2 + I 2 Q 1 ) 


parameters (tp,tq) are increased. For the TKO applied to the 
signal, at an input SNR of 17 dB, the output SNRs for the 
operators thg, thg and respectively (a) 6.65 dB, (b) 

24.92 dB, and (c) 29.52 dB. Evidently, the four-sample input 
TKO provides a better noise performance compared to a 
three-sample input TKO Ihg. This can be attributed to the fact 
that the mean and variance of the noise terms are larger in the 
latter case. Similar remarks hold for the TKO applied to the 
signal derivative. 

Furthermore, depending upon the structure of the operator 
kernel J and the input noise correlation, a finite number of 
eigenvalues associated with the product M J in (5) are always 
found to be negative. For instance, at an input SNR of 17 dB, 
the eigenvalues of MJ computed for the operators T'g, ipQ and 
4^4 are respectively, (.836,.475,-.361), (.952,.498,-.455), 
and (.655,.376,—.321,—.042). Negative eigenvalues corre¬ 
spond physically to the fact that the operator output V(t) can 
go negative for some time f (> 0), although there can still be 
a positive mean value. If the delay parameters (tp,tq) are large 
(vis-a-vis the correlation time of the input noise, 6c), we have 
a large number of eigenvalues and an asymptotically normal 
density function can be expected. Conversely, for small delay 
parameters we expect large departures from normality, as is 
evidenced in the plots in Fig. 2. From Fig. 2, it is also evident 
that the probability of excursion towards the negative region 
is lower for T'l compared to Ihg and Ikj. 

For narrowband inputs, ratios of quadratic forms that yield 
the sine and tangent of the phase difference are respec¬ 
tively, (a) 2 {IiQ 2 - hQi)l{Ii + I2 + Qi + Ql) (b) 
{I1Q2 - l2Qi)/iIiQ2 + hQl), where h and Qi are the 
in-phase and quadrature components, i = 1,2 denoting two 
adjacent samples. Fig. 3 shows the density functions of the 
individual quadratic forms and their ratios. For the quadratic 
form {I1Q2 — hQl), the mean and the higher odd order 
moments of the noise term are zero. It is easy to show that the 
covariance of {IiQ 2 — hQi) and {I 1 +I 2 +Q 1 +Q 2 ) is greater 
than that of {I1Q2 — hQi) and {I1Q2 + hQi)- This gives a 
marginally higher variance for the ratio in (b) compared to that 
in (a). Pdf of the ratios shown in the inset of Fig. 3 confirms 
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Fig. 4. TKO output before and after post-detection filtering using a three- 
point Binomial filter with an impulse response (1,2,1). Input SNR =11 dB. 


these observations. 

A way of increasing the signal mean value and ensuring 
the positivity of the operator outputs is to sum several TKO 
outputs. If the observation time is long enough to allow for 
decorrelation in time, the constituent outputs are independent 
and the cumulants add. This ensures an increase in the 
mean value. Positivity can also be achieved by introducing 
a spectrally narrow low-pass filter at the TKO post-detection 
stage (e.g., a filter with a long time constant compared to 
6c)- Integration attenuates high frequency spectral components 
and is an effective option for filtering the strong time-varying 
fluctuations induced by the signalx noise terms. A long time 
constant however compromises the “instantaneous” property of 
the TKO. Results of applying a three-point binomial filter at 
the TKO output 'I'[i(f)] is shown in Fig. 4. The short impulse 
response (1,2,1) of the filter preserves the superior localization 
property of the TKO and restores positivity of the operator 
outputs. When filtering is not an option, simple thresholding 
to remove negative outputs or outputs near zero is a viable 
alternative that can retain the speed advantage of the operator, 
if the penalty in terms of missed data is not a concern. 
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